TRAF3 loss-of-function reveals the noncanonical NF-κB pathway as a therapeutic target in diffuse large B cell lymphoma.

Here, we report recurrent focal deletions of the chr14q32.31-32 locus, including TRAF3, a negative regulator of NF-κB signaling, in de novo diffuse large B cell lymphoma (DLBCL) (24/324 cases). Integrative analysis revealed an association between TRAF3 copy number loss with accumulation of NIK, the central noncanonical (NC) NF-κB kinase, and increased NC NF-κB pathway activity. Accordingly, TRAF3 genetic ablation in isogenic DLBCL model systems caused upregulation of NIK and enhanced NC NF-κB downstream signaling. Knockdown or pharmacological inhibition of NIK in TRAF3-deficient cells differentially impaired their proliferation and survival, suggesting an acquired onco-addiction to NC NF-κB. TRAF3 ablation also led to exacerbated secretion of the immunosuppressive cytokine IL-10. Coculturing of TRAF3-deficient DLBCL cells with CD8+ T cells impaired the induction of Granzyme B and interferon (IFN) γ, which were restored following neutralization of IL-10. Our findings corroborate a direct relationship between TRAF3 genetic alterations and NC NF-κB activation, and highlight NIK as a potential therapeutic target in a defined subset of DLBCL.

Spatially Resolved Tumor Microenvironment Predicts Treatment Outcomes in Relapsed/Refractory Hodgkin Lymphoma.

PURPOSE: About a third of patients with relapsed or refractory classic Hodgkin lymphoma (r/r CHL) succumb to their disease after high-dose chemotherapy followed by autologous stem-cell transplantation (HDC/ASCT). Here, we aimed to describe spatially resolved tumor microenvironment (TME) ecosystems to establish novel biomarkers associated with treatment failure in r/r CHL. PATIENTS AND METHODS: We performed imaging mass cytometry (IMC) on 71 paired primary diagnostic and relapse biopsies using a marker panel specific to CHL biology. For each cell type in the TME, we calculated a spatial score measuring the distance of nearest neighbor cells to the malignant Hodgkin Reed Sternberg cells within the close interaction range. Spatial scores were used as features in prognostic model development for post-ASCT outcomes. RESULTS: Highly multiplexed IMC data revealed shared TME patterns in paired diagnostic and early r/r CHL samples, whereas TME patterns were more divergent in pairs of diagnostic and late relapse samples. Integrated analysis of IMC and single-cell RNA sequencing data identified unique architecture defined by CXCR5+ Hodgkin and Reed Sternberg (HRS) cells and their strong spatial relationship with CXCL13+ macrophages in the TME. We developed a prognostic assay (RHL4S) using four spatially resolved parameters, CXCR5+ HRS cells, PD1+CD4+ T cells, CD68+ tumor-associated macrophages, and CXCR5+ B cells, which effectively separated patients into high-risk versus low-risk groups with significantly different post-ASCT outcomes. The RHL4S assay was validated in an independent r/r CHL cohort using a multicolor immunofluorescence assay. CONCLUSION: We identified the interaction of CXCR5+ HRS cells with ligand-expressing CXCL13+ macrophages as a prominent crosstalk axis in relapsed CHL. Harnessing this TME biology, we developed a novel prognostic model applicable to r/r CHL biopsies, RHL4S, opening new avenues for spatial biomarker development.

Distinct Gene Expression Patterns of Calcium Channels and Related Signaling Pathways Discovered in Lymphomas.

Cell surface calcium (Ca2+) channels permit Ca2+ ion influx, with Ca2+ taking part in cellular functions such as proliferation, survival, and activation. The expression of voltage-dependent Ca2+ (CaV) channels may modulate the growth of hematologic cancers. Profile analysis of Ca2+ channels, with a focus on the Ca2+ release-activated Ca2+ (CRAC) and L-type CaV channels, was performed on RNA sequencing data from lymphoma cell lines and samples derived from patients with diffuse large B cell lymphoma (DLBCL). CaV1.2 expression was found to be elevated in classical Hodgkin lymphoma (CHL) cell lines when compared to other B cell lymphoma cell lines. In contrast, CHL exhibited reduced expression of ORAI2 and STIM2. In our differential expression analysis comparing activated B cell-like DLBCL (ABC-DLBCL) and germinal centre B cell-like DLBCL (GCB-DLBCL) patient samples, ABC-DLBCL revealed stronger expression of CaV1.3, whereas CaV1.1, CaV1.2, and CaV1.4 showed greater expression levels in GCB-DLBCL. Interestingly, no differences in ORAI/STIM expression were noted in the patient samples. As Ca2+ is known to bind to calmodulin, leading to calcineurin activation and the passage of nuclear factor of activated T cells (NFAT) to the cell nucleus, pathways for calcineurin, calmodulin, NFAT, and Ca2+ signaling were also analyzed by gene set enrichment analysis. The NFAT and Ca2+ signaling pathways were found to be upregulated in the CHL cell lines relative to other B cell lymphoma cell lines. Furthermore, the calmodulin and Ca2+ signaling pathways were shown to be downregulated in the ABC-DLBCL patient samples. The findings of this study suggest that L-type CaV channels and Ca2+-related pathways could serve as differentiating components for biologic therapies in targeted lymphoma treatments.

Tumor-associated antigen PRAME exhibits dualistic functions that are targetable in diffuse large B cell lymphoma.

PRAME is a prominent member of the cancer testis antigen family of proteins, which triggers autologous T cell-mediated immune responses. Integrative genomic analysis in diffuse large B cell lymphoma (DLBCL) uncovered recurrent and highly focal deletions of 22q11.22, including the PRAME gene, which were associated with poor outcome. PRAME-deleted tumors showed cytotoxic T cell immune escape and were associated with cold tumor microenvironments. In addition, PRAME downmodulation was strongly associated with somatic EZH2 Y641 mutations in DLBCL. In turn, PRC2-regulated genes were repressed in isogenic PRAME-KO lymphoma cell lines, and PRAME was found to directly interact with EZH2 as a negative regulator. EZH2 inhibition with EPZ-6438 abrogated these extrinsic and intrinsic effects, leading to PRAME expression and microenvironment restoration in vivo. Our data highlight multiple functions of PRAME during lymphomagenesis and provide a preclinical rationale for synergistic therapies combining epigenetic reprogramming with PRAME-targeted therapies.

Single-cell profiling reveals the importance of CXCL13/CXCR5 axis biology in lymphocyte-rich classic Hodgkin lymphoma.

Lymphocyte-rich classic Hodgkin lymphoma (LR-CHL) is a rare subtype of Hodgkin lymphoma. Recent technical advances have allowed for the characterization of specific cross-talk mechanisms between malignant Hodgkin Reed-Sternberg (HRS) cells and different normal immune cells in the tumor microenvironment (TME) of CHL. However, the TME of LR-CHL has not yet been characterized at single-cell resolution. Here, using single-cell RNA sequencing (scRNA-seq), we examined the immune cell profile of 8 cell suspension samples of LR-CHL in comparison to 20 samples of the mixed cellularity (MC, 9 cases) and nodular sclerosis (NS, 11 cases) subtypes of CHL, as well as 5 reactive lymph node controls. We also performed multicolor immunofluorescence (MC-IF) on tissue microarrays from the same patients and an independent validation cohort of 31 pretreatment LR-CHL samples. ScRNA-seq analysis identified a unique CD4+ helper T cell subset in LR-CHL characterized by high expression of Chemokine C-X-C motif ligand 13 (CXCL13) and PD-1. PD-1+CXCL13+ T cells were significantly enriched in LR-CHL compared to other CHL subtypes, and spatial analyses revealed that in 46% of the LR-CHL cases these cells formed rosettes surrounding HRS cells. MC-IF analysis revealed CXCR5+ normal B cells in close proximity to CXCL13+ T cells at significantly higher levels in LR-CHL. Moreover, the abundance of PD-1+CXCL13+ T cells in the TME was significantly associated with shorter progression-free survival in LR-CHL (P = 0.032). Taken together, our findings strongly suggest the pathogenic importance of the CXCL13/CXCR5 axis and PD-1+CXCL13+ T cells as a treatment target in LR-CHL.

Combined EZH2 Inhibition and IKAROS Degradation Leads to Enhanced Antitumor Activity in Diffuse Large B-cell Lymphoma.

PURPOSE: The efficacy of EZH2 inhibition has been modest in the initial clinical exploration of diffuse large B-cell lymphoma (DLBCL), yet EZH2 inhibitors are well tolerated. Herein, we aimed to uncover genetic and pharmacologic opportunities to enhance the clinical efficacy of EZH2 inhibitors in DLBCL. EXPERIMENTAL DESIGN: We conducted a genome-wide sensitizing CRISPR/Cas9 screen with tazemetostat, a catalytic inhibitor of EZH2. The sensitizing effect of IKZF1 loss of function was then validated and leveraged for combination treatment with lenalidomide. RNA sequencing (RNA-seq) and chromatin immunoprecipitation sequencing analyses were performed to elucidate transcriptomic and epigenetic changes underlying synergy. RESULTS: We identified IKZF1 knockout as the top candidate for sensitizing DLBCL cells to tazemetostat. Treating cells with tazemetostat and lenalidomide, an immunomodulatory drug that selectively degrades IKAROS and AIOLOS, phenocopied the effects of the CRISPR/Cas9 screen. The combined drug treatment triggered either cell-cycle arrest or apoptosis in a broad range of DLBCL cell lines, regardless of EZH2 mutational status. Cell-line-based xenografts also showed slower tumor growth and prolonged survival in the combination treatment group. RNA-seq analysis revealed strong upregulation of interferon signaling and antiviral immune response signatures. Gene expression of key immune response factors such as IRF7 and DDX58 were induced in cells treated with lenalidomide and tazemetostat, with a concomitant increase of H3K27 acetylation at their promoters. Furthermore, transcriptome analysis demonstrated derepression of endogenous retroviruses after combination treatment. CONCLUSIONS: Our data underscore the synergistic interplay between IKAROS degradation and EZH2 inhibition on modulating epigenetic changes and ultimately enhancing antitumor effects in DLBCL.

Single-cell Profiles and Prognostic Impact of Tumor-Infiltrating Lymphocytes Coexpressing CD39, CD103, and PD-1 in Ovarian Cancer.

PURPOSE: Tumor-infiltrating lymphocytes (TIL) are strongly associated with survival in most cancers; however, the tumor-reactive subset that drives this prognostic effect remains poorly defined. CD39, CD103, and PD-1 have been independently proposed as markers of tumor-reactive CD8+ TIL in various cancers. We evaluated the phenotype, clonality, and prognostic significance of TIL expressing various combinations of these markers in high-grade serous ovarian cancer (HGSC), a malignancy in need of more effective immunotherapeutic approaches. EXPERIMENTAL DESIGN: Expression of CD39, CD103, PD-1, and other immune markers was assessed by high-dimensional flow cytometry, single-cell sequencing, and multiplex immunofluorescence of primary and matched pre/post-chemotherapy HGSC specimens. RESULTS: Coexpression of CD39, CD103, and PD-1 ("triple-positive" phenotype) demarcated subsets of CD8+ TIL and CD4+ regulatory T cells (Treg) with a highly activated/exhausted phenotype. Triple-positive CD8+ TIL exhibited reduced T-cell receptor (TCR) diversity and expressed genes involved in both cytolytic and humoral immunity. Triple-positive Tregs exhibited higher TCR diversity and a tumor-resident phenotype. Triple-positive TIL showed superior prognostic impact relative to TIL expressing other combinations of these markers. TIGIT was uniquely upregulated on triple-positive CD8+ effector cells relative to their CD4+ Treg counterparts. CONCLUSIONS: Coexpression of CD39, CD103, and PD-1 demarcates highly activated CD8+ and CD4+ TIL with inferred roles in cytolytic, humoral, and regulatory immune functions. Triple-positive TIL demonstrate exceptional prognostic significance and express compelling targets for combination immunotherapy, including PD-1, CD39, and TIGIT.

Mutational landscape of gray zone lymphoma.

The mutational landscape of gray zone lymphoma (GZL) has not yet been established, and differences from related entities are largely unknown. Here, we studied coding sequence mutations of 50 Epstein-Barr virus (EBV)-negative GZLs and 20 polymorphic EBV+ diffuse large B-cell lymphoma (DLBCL) not otherwise specified (poly-EBV-L) in comparison with classical Hodgkin lymphoma (cHL), primary mediastinal large B-cell lymphoma (PMBCL), and DLBCL. Exomes of 21 GZL and 7 poly-EBV-L cases, along with paired constitutional DNA, were analyzed as a discovery cohort, followed by targeted sequencing of 217 genes in an extension cohort of 29 GZL and 13 poly-EBV-L cases. GZL cases with thymic niche involvement (anterior mediastinal mass) exhibited a mutation profile closely resembling cHL and PMBCL, with SOCS1 (45%), B2M (45%), TNFAIP3 (35%), GNA13 (35%), LRRN3 (32%), and NFKBIA (29%) being the most recurrently mutated genes. In contrast, GZL cases without thymic niche involvement (n = 18) had a significantly distinct pattern that was enriched in mutations related to apoptosis defects (TP53 [39%], BCL2 [28%], BIRC6 [22%]) and depleted in GNA13, XPO1, or NF-κB signaling pathway mutations (TNFAIP3, NFKBIE, IKBKB, NFKBIA). They also exhibited more BCL2/BCL6 rearrangements compared with thymic GZL. Poly-EBV-L cases presented a distinct mutational profile, including STAT3 mutations and a significantly lower coding mutation load in comparison with EBV- GZL. Our study highlights characteristic mutational patterns in GZL associated with presentation in the thymic niche, suggesting a common cell of origin and disease evolution overlapping with related anterior mediastinal lymphomas.

Single-Cell Transcriptome Analysis Reveals Disease-Defining T-cell Subsets in the Tumor Microenvironment of Classic Hodgkin Lymphoma.

Hodgkin lymphoma is characterized by an extensively dominant tumor microenvironment (TME) composed of different types of noncancerous immune cells with rare malignant cells. Characterization of the cellular components and their spatial relationship is crucial to understanding cross-talk and therapeutic targeting in the TME. We performed single-cell RNA sequencing of more than 127,000 cells from 22 Hodgkin lymphoma tissue specimens and 5 reactive lymph nodes, profiling for the first time the phenotype of the Hodgkin lymphoma-specific immune microenvironment at single-cell resolution. Single-cell expression profiling identified a novel Hodgkin lymphoma-associated subset of T cells with prominent expression of the inhibitory receptor LAG3, and functional analyses established this LAG3+ T-cell population as a mediator of immunosuppression. Multiplexed spatial assessment of immune cells in the microenvironment also revealed increased LAG3+ T cells in the direct vicinity of MHC class II-deficient tumor cells. Our findings provide novel insights into TME biology and suggest new approaches to immune-checkpoint targeting in Hodgkin lymphoma. SIGNIFICANCE: We provide detailed functional and spatial characteristics of immune cells in classic Hodgkin lymphoma at single-cell resolution. Specifically, we identified a regulatory T-cell-like immunosuppressive subset of LAG3+ T cells contributing to the immune-escape phenotype. Our insights aid in the development of novel biomarkers and combination treatment strategies targeting immune checkpoints.See related commentary by Fisher and Oh, p. 342.This article is highlighted in the In This Issue feature, p. 327.

Integrative genomic analysis identifies key pathogenic mechanisms in primary mediastinal large B-cell lymphoma.

Primary mediastinal large B-cell lymphoma (PMBL) represents a clinically and pathologically distinct subtype of large B-cell lymphomas. Furthermore, molecular studies, including global gene expression profiling, have provided evidence that PMBL is more closely related to classical Hodgkin lymphoma (cHL). Although targeted sequencing studies have revealed a number of mutations involved in PMBL pathogenesis, a comprehensive description of disease-associated genetic alterations and perturbed pathways is still lacking. Here, we performed whole-exome sequencing of 95 PMBL tumors to inform on oncogenic driver genes and recurrent copy number alterations. The integration of somatic gene mutations with gene expression signatures provides further insights into genotype-phenotype interrelation in PMBL. We identified highly recurrent oncogenic mutations in the Janus kinase-signal transducer and activator of transcription and nuclear factor κB pathways, and provide additional evidence of the importance of immune evasion in PMBL (CIITA, CD58, B2M, CD274, and PDCD1LG2). Our analyses highlight the interferon response factor (IRF) pathway as a putative novel hallmark with frequent alterations in multiple pathway members (IRF2BP2, IRF4, and IRF8). In addition, our integrative analysis illustrates the importance of JAK1, RELB, and EP300 mutations driving oncogenic signaling. The identified driver genes were significantly more frequently mutated in PMBL compared with diffuse large B-cell lymphoma, whereas only a limited number of genes were significantly different between PMBL and cHL, emphasizing the close relation between these entities. Our study, performed on a large cohort of PMBL, highlights the importance of distinctive genetic alterations for disease taxonomy with relevance for diagnostic evaluation and therapeutic decision-making.

BPG: Seamless, automated and interactive visualization of scientific data.

BACKGROUND: We introduce BPG, a framework for generating publication-quality, highly-customizable plots in the R statistical environment. RESULTS: This open-source package includes multiple methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it suitable for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for integration with computational pipelines. CONCLUSION: BPG provides a new approach for linking interactive and scripted data visualization and is available at http://labs.oicr.on.ca/boutros-lab/software/bpg or via CRAN at https://cran.r-project.org/web/packages/BoutrosLab.plotting.general.

High-resolution architecture and partner genes of MYC rearrangements in lymphoma with DLBCL morphology.

Genomic rearrangements in the MYC locus occur in ∼12% of lymphomas with diffuse large B-cell lymphoma (DLBCL) morphology and are associated with inferior outcome. Previous studies exploring MYC rearrangements have primarily used fluorescence in situ hybridization (FISH) assays to characterize break-apart status but have rarely examined breakpoint location, and in some cases have not examined partner identity. We performed targeted sequencing of MYC, BCL2, BCL6, and the immunoglobulin (IG) loci in 112 tumors with DLBCL morphology harboring MYC rearrangement. We characterized the location of the MYC rearrangement at base pair resolution and identified the partner in 88 cases. We observed a cluster of breakpoints upstream of the MYC coding region and in intron 1 (the "genic cluster"). Genic cluster rearrangements were enriched for translocations involving IGH (80%), whereas nongenic rearrangements occurred mostly downstream of the MYC gene with a variety of partners, including IGL and IGK Other recurrent partners included BCL6, ZCCHC7, and RFTN1, which has not previously been described as a MYC partner. We compared 2 commercially available FISH break-apart assays for the MYC locus and observed discordant results in 32% of cases examined, including some with MYC-IGL and MYC-IGK rearrangements. In cases of high-grade B-cell lymphoma with MYC and BCL2 and/or BCL6 rearrangement (HGBL-DH), so-called "double-hit" lymphomas, the majority of MYC rearrangements had non-IG partners (65%), with breakpoints outside the genic cluster (72%). In patients with de novo HGBL-DH of DLBCL morphology, MYC-IG rearrangements showed a trend toward inferior time to progression and overall survival compared with MYC-non-IG rearrangements. Our data reveal clinically relevant architecture of MYC rearrangements in lymphomas with DLBCL morphology.

Compendium of TCDD-mediated transcriptomic response datasets in mammalian model systems.

BACKGROUND: 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent congener of the dioxin class of environmental contaminants. Exposure to TCDD causes a wide range of toxic outcomes, ranging from chloracne to acute lethality. The severity of toxicity is highly dependent on the aryl hydrocarbon receptor (AHR). Binding of TCDD to the AHR leads to changes in transcription of numerous genes. Studies evaluating the transcriptional changes brought on by TCDD may provide valuable insight into the role of the AHR in human health and disease. We therefore compiled a collection of transcriptomic datasets that can be used to aid the scientific community in better understanding the transcriptional effects of ligand-activated AHR. RESULTS: Specifically, we have created a datasets package - TCDD.Transcriptomics - for the R statistical environment, consisting of 63 unique experiments comprising 377 samples, including various combinations of 3 species (human derived cell lines, mouse and rat), 4 tissue types (liver, kidney, white adipose tissue and hypothalamus) and a wide range of TCDD exposure times and doses. These datasets have been fully standardized using consistent preprocessing and annotation packages (available as of September 14, 2015). To demonstrate the utility of this R package, a subset of "AHR-core" genes were evaluated across the included datasets. Ahrr, Nqo1 and members of the Cyp family were significantly induced following exposure to TCDD across the studies as expected while Aldh3a1 was induced specifically in rat liver. Inmt was altered only in liver tissue and primarily by rat-AHR. CONCLUSIONS: Analysis of the "AHR-core" genes demonstrates a continued need for studies surrounding the impact of AHR-activity on the transcriptome; genes believed to be consistently regulated by ligand-activated AHR show surprisingly little overlap across species and tissues. Until now, a comprehensive assessment of the transcriptome across these studies was challenging due to differences in array platforms, processing methods and annotation versions. We believe that this package, which is freely available for download ( http://labs.oicr.on.ca/boutros-lab/tcdd-transcriptomics ) will prove to be a highly beneficial resource to the scientific community evaluating the effects of TCDD exposure as well as the variety of functions of the AHR.

Comprehensive characterization of programmed death ligand structural rearrangements in B-cell non-Hodgkin lymphomas.

Programmed death ligands (PDLs) are immune-regulatory molecules that are frequently affected by chromosomal alterations in B-cell lymphomas. Although PDL copy-number variations are well characterized, a detailed and comprehensive analysis of structural rearrangements (SRs) and associated phenotypic consequences is largely lacking. Here, we used oligonucleotide capture sequencing of 67 formalin-fixed paraffin-embedded tissues derived from primary B-cell lymphomas and 1 cell line to detect and characterize, at base-pair resolution, SRs of the PDL locus (9p24.1; harboring PDL1/CD274 and PDL2/PDCD1LG2). We describe 36 novel PDL SRs, including 17 intrachromosomal events (inversions, duplications, deletions) and 19 translocations involving BZRAP-AS1, CD44, GET4, IL4R, KIAA0226L, MID1, RCC1, PTPN1 and segments of the immunoglobulin loci. Moreover, analysis of the precise chromosomal breakpoints reveals 2 distinct cluster breakpoint regions (CBRs) within either CD274 (CBR1) or PDCD1LG2 (CBR2). To determine the phenotypic consequences of these SRs, we performed immunohistochemistry for CD274 and PDCD1LG2 on primary pretreatment biopsies and found that PDL SRs are significantly associated with PDL protein expression. Finally, stable ectopic expression of wild-type PDCD1LG2 and the PDCD1LG2-IGHV7-81 fusion showed, in coculture, significantly reduced T-cell activation. Taken together, our data demonstrate the complementary utility of fluorescence in situ hybridization and capture sequencing approaches and provide a classification scheme for PDL SRs with implications for future studies using PDL immune-checkpoint inhibitors in B-cell lymphomas.

Spatial genomic heterogeneity within localized, multifocal prostate cancer.

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

Cross-species transcriptomic analysis elucidates constitutive aryl hydrocarbon receptor activity.

BACKGROUND: Research on the aryl hydrocarbon receptor (AHR) has largely focused on variations in toxic outcomes resulting from its activation by halogenated aromatic hydrocarbons. But the AHR also plays key roles in regulating pathways critical for development, and after decades of research the mechanisms underlying physiological regulation by the AHR remain poorly characterized. Previous studies identified several core genes that respond to xenobiotic AHR ligands across a broad range of species and tissues. However, only limited inferences have been made regarding its role in regulating constitutive gene activity, i.e. in the absence of exogenous ligands. To address this, we profiled transcriptomic variations between AHR-active and AHR-less-active animals in the absence of an exogenous agonist across five tissues, three of which came from rats (hypothalamus, white adipose and liver) and two of which came from mice (kidney and liver). Because AHR status alone has been shown sufficient to alter transcriptomic responses, we reason that by contrasting profiles amongst AHR-variant animals, we may elucidate effects of the AHR on constitutive mRNA abundances. RESULTS: We found significantly more overlap in constitutive mRNA abundances amongst tissues within the same species than from tissues between species and identified 13 genes (Agt, Car3, Creg1, Ctsc, E2f6, Enpp1, Gatm, Gstm4, Kcnj8, Me1, Pdk1, Slc35a3, and Sqrdl) that are affected by AHR-status in four of five tissues. One gene, Creg1, was significantly up-regulated in all AHR-less-active animals. We also find greater overlap between tissues at the pathway level than at the gene level, suggesting coherency to the AHR signalling response within these processes. Analysis of regulatory motifs suggests that the AHR mostly mediates transcriptional regulation via direct binding to response elements. CONCLUSIONS: These findings, though preliminary, present a platform for further evaluating the role of the AHR in regulation of constitutive mRNA levels and physiologic function.

Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study.

BACKGROUND: Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. METHODS: We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. FINDINGS: Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. INTERPRETATION: This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. FUNDING: Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

SeqControl: process control for DNA sequencing.

As high-throughput sequencing continues to increase in speed and throughput, routine clinical and industrial application draws closer. These 'production' settings will require enhanced quality monitoring and quality control to optimize output and reduce costs. We developed SeqControl, a framework for predicting sequencing quality and coverage using a set of 15 metrics describing overall coverage, coverage distribution, basewise coverage and basewise quality. Using whole-genome sequences of 27 prostate cancers and 26 normal references, we derived multivariate models that predict sequencing quality and depth. SeqControl robustly predicted how much sequencing was required to reach a given coverage depth (area under the curve (AUC) = 0.993), accurately classified clinically relevant formalin-fixed, paraffin-embedded samples, and made predictions from as little as one-eighth of a sequencing lane (AUC = 0.967). These techniques can be immediately incorporated into existing sequencing pipelines to monitor data quality in real time. SeqControl is available at http://labs.oicr.on.ca/Boutros-lab/software/SeqControl/.

Validating reference genes within a mouse model system of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) toxicity.

BACKGROUND: Quantitative real-time PCR (qPCR) is the "gold-standard" technique for measuring mRNA abundances. To correctly compare samples and generate biologically valid results, qPCR data usually require comprehensive normalization to account for sample content variation between reactions. The most common normalization approaches use one or more endogenous controls (reference or house-keeping genes) to adjust the measured levels of experimental genes appropriately. Ideal reference genes are those that display minimal variation across experimental conditions, and thus can vary widely across different biological systems. In particular, toxicogenomic studies of transcriptionally-disruptive toxins, like 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), require careful consideration of reference genes. RESULTS: We examined seven candidate reference genes in 199 mice varying in genotype and time/dose of TCDD exposure. We assessed gene-stability in four ways: (1) the variance of the raw Cq values across biological replicates, (2) the fold-change from basal mRNA levels following treatment, (3) the inter- and intra-group stability evaluated using the NormFinder algorithm, (4) the comparative ΔCq method for each candidate gene. Univariate analyses showed Hprt and Eef1a1 are the two most stable individual reference genes. It has been suggested that using multiple genes would produce a more consistent normalization factor; multivariate analysis was performed using NormFinder. In general, stability increased with the number of genes used, but specific gene-combinations synergized. CONCLUSIONS: We have validated seven reference genes for use in analyzing mRNA abundances in mouse models of TCDD toxicity. The use of multiple reference genes increases stability, providing more consistent normalization and more reliable results. The number of reference genes used should be maximized, based on experimental capabilities (platform, sample availability, etc.). Our results show the benefit of validating reference genes using multiple methods prior to generating large biological datasets.

Inter-strain heterogeneity in rat hepatic transcriptomic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD).

The biochemical and toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) have been the subject of intense study for decades. It is now clear that essentially all TCDD-induced toxicities are mediated by DNA-protein interactions involving the Aryl Hydrocarbon Receptor (AHR). Nevertheless, it remains unknown which AHR target genes cause TCDD toxicities. Several groups, including our own, have developed rodent model systems to probe these questions. mRNA expression profiling of these model systems has revealed significant inter-species heterogeneity in rodent hepatic responses to TCDD. It has remained unclear if this variability also exists within a species, amongst rodent strains. To resolve this question, we profiled the hepatic transcriptomic response to TCDD of diverse rat strains (L-E, H/W, F344 and Wistar rats) and two lines derived from L-E×H/W crosses, at consistent age, sex, and dosing (100 µg/kg TCDD for 19 h). Using this uniquely consistent dataset, we show that the majority of TCDD-induced alterations in mRNA abundance are strain/line-specific: only 11 genes were affected by TCDD across all strains, including well-known dioxin-responsive genes such as Cyp1a1 and Nqo1. Our analysis identified two novel universally dioxin-responsive genes as well as 4 genes induced by TCDD in dioxin-sensitive rats only. These 6 genes are strong candidates to explain TCDD-related toxicities, so we validated them using 152 animals in time-course (0 to 384 h) and dose-response (0 to 3000 µg/kg) experiments. This study reveals that different rat strains exhibit dramatic transcriptional heterogeneity in their hepatic responses to TCDD and that inter-strain comparisons can help identify candidate toxicity-related genes.